Field induced transitions in a kagome antiferromagnet 
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The thermal order by disorder effect in magnetic field is studied for a classical Heisenberg antiferro- 
magnet on the kagome lattice. Using analytical arguments we predict a unique H—T phase diagram 
for this strongly frustrated magnet: states with a coplanar and a uniaxial triatic order parameters 
respectively at low and high magnetic fields and an incompressible collinear spin-liquid state at a 
one-third of the saturation field. We also present the Monte Carlo data which confirm existence of 
these phases. 
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Geometrical frustration in lattice spin models is re- 
sponsible for complete suppression of conventional mag- 
netic order and appearance of nonmagnetic spin-liquids 
or states with exotic order parameters |lj . Applied mag- 
netic field can further enhance frustration. A well known 
example is a weakly frustrated triangular lattice antifer- 
romagnet, which acquires an additional continuous de- 
generacy in external field g] . Investigation of high- field 
effects in strongly frustrated magnets poses a new chal- 
lenge for experimental and theoretical studies pL|[ . 

The Heisenberg antiferromagnet on the kagome lattice 
(see Fig. [l]) is a strongly frustrated spin model, which is 
approximately realized in a number of insulating layered 
magnets: SrCr 8 Ga40ig || and Ba2Sn 2 Ga3ZnCr7022 Jq] 
(both with S = 3/2), KFe 3 (OH) 6 (S0 4 )2 (S = 5/2) [§. 
Gd3Ga50i2 (S — 7/2) is another frustrated magnet on a 
related three-dimensional garnet lattice of corner-sharing 
triangles, which is often called a hyper-kagome lattice. 
This magnet has a weak exchange constant J~ IK and 
a peculiar unexplained phase diagram in magnetic field 
H . Motivated by the above materials with large values of 
spin we investigate in this Letter the finite-temperature 
magnetization process of a classical antiferromagnet on 
the kagome lattice. We predict three distinctive field 
regimes below the saturation field i? sa t, where exotic spin 
phases are stabilized. 

The Hamiltonian of a nearest-neighbor Heisenberg an- 
tiferromagnet on the kagome lattice with classical unit 
spins can be written up to a constant term as: 






(1) 



where the sum runs over all triangles, Sa is the total spin 
of a triangular plaquette and N& = | N is the number of 
plaquettes on an iV-site lattice. The zero-field classical 
constraint Sa = fails to define a unique ground state. 
The ground state coplanar configurations are constructed 
by fixing spins on a first plaquette to a 120°-structure 
with left or right chirality for the triad e a , e&, and e c 
and, then, tiling this triad over the whole lattice in a way 
that three spins on every plaquette are different [|| fLUf . 



The number of all such states for the kagome lattice is 
known exactly: U N with U = 1.1347... |ll]]. Nonplanar 
ground states are constructed from planar configurations 
by identifying so-called weathervane defects |§ . 

The thermal order by disorder effect can appear be- 
cause of a different entropy of short-wavelength fluctua- 
tions above degenerate configurations |L2J. All coplanar 
states for a kagome antiferromagnet show identical har- 
monic spectra. They have one flat zero-energy branch 
with N4 = hN modes, which corresponds to anharmonic 
quartic excitations flsHlCfl . A nonzero harmonic mode de- 
scribed by a classical coordinate y has energy increase 
AE 2 ~ Jy 2 and contributes \T ln( J/T) to the ther- 
modynamic potential, whereas a soft quartic mode with 
AE4 ~ Jy 4 makes a reduced contribution of jTln(J/T). 
Coplanar configurations have the largest number of soft 
modes and are, therefore, selected by thermal fluctua- 
tions H . Soft modes in coplanar states correspond to al- 
ternate tilting of spins out of the ground state spin plane 
around elementary hexagons. There are ^N hexagons 
on the kagome lattice. The counting of soft modes from 
such a geometrical point of view, thus, agrees with the 
spin- wave analysis. 

Harmonic fluctuations do not select between vari- 
ous planar configurations, though the highest statistical 
weight corresponds to a so-called \/3 x %/3 structure [gjj . 
The low-temperature phase has nematic correlations of 
the chirality vectors defined on every plaquette ||. Be- 
sides, the broken symmetries include selection of the spin 
triad. The residual symmetry group is determined by the 
two elements: CqT and U2, Cg is a spin rotation by 7r/3 
about normal to the spin plane, T is the time-reversal, 
and U2 is a rotation by 7r about one of the three in-plane 
axes. The order parameter for low-temperature phases of 
the kagome antiferromagnet is a three-point correlation 
function taken on one plaquette: 



(S a {t 1 )S^(T 2 )S^(T i ))=S a ^ 



(2) 



which has only a spin part described by a fully symmetric 
traceless tensor S a ^ P,M ■ For coplanar states at if = 
this tensor is parameterized as 



Qa^-y 



&n Ch&r 



(3) 



Note, that the three-spin (triatic) order parameter breaks 
explicitly the time-reversal symmetry, which determines 
its nontrivial coupling to an applied field. 

In a finite magnetic field the classical energy Eq. (fil) 
reaches the minimum provided that 



'A 



H/(2J) 



(4) 



for every plaquette. Since |Sa| < 3, above the saturation 
field H sat — 6 J all spins are aligned parallel to H. We 
are interested in the field range < H < H sat , where the 
classical ground state possesses an infinite degeneracy. 
The degeneracy is further enhanced by a field because 
the minimum energy constraint does not require spins 
on one plaquette to lie in the same plane. Shcnder and 
Holdsworth made the only attempt to understand the 
finite field behavior of the kagome antiferromagnet [[Wj . 
They noticed that spin coplanarity responsible for the 
increased number of soft modes is preserved if the spin 
plane is parallel to an external field. Hence, thermal 
fluctuations create magnetic anisotropy, which orients at 
low temperatures the spin plane. 

There is an additional effect missed in the above con- 
sideration and related to the triatic nature of the spin 
order parameter. Once selection of the spin plane takes 
place, an extra macroscopic degree of freedom appears: 
orientation of the spin triad inside the plane relative to 
H. At T = all orientations have the same classical 
energy, though the spin triad distorts differently, when 
it forms different angles p with the magnetic field. This 
distortion modifies harmonic fluctuations inside the spin 
plane. In such cases with a one-parameter continuous de- 
generacy thermal fluctuations always produce the order 
by disorder effect and select a homogeneous state for one 
particular value of the parameter |Q . The orientational 
field effect for the triatic order parameter is described by 
an invariant £ an ~ H a H fj ' H~< S a ^ ex H 3 cos Sip, which 
is present in the Landau functional because of a broken 
time-reversal symmetry in the triatic phase. The mini- 
mum energy is reached for p = or -k depending on a 
sign of the prefactor, i.e., one of the triad vectors has to 
be either parallel or antiparallel to the applied field. 

To find the equilibrium orientation of the spin triad 
in external field, we performed spin- wave calculations. 
The harmonic spectrum of a coplanar state at H ^ 
depends not only on the triad orientation but also on a 
way the triad is tiled over the lattice. We used, there- 
fore, the V3 x \/3 and the q — structures (l(| as two 
example cases. For both structures p — ir is selected 
by fluctuations with E an = 5 .0 T (H / H sat ) 3 N cosSp and 
1.3 T(H/H sat ) 3 N cos 3p, respectively. Thus, thermal 
fluctuations in a weak magnetic field create anisotropy 
for the spin plane as well as inside the plane. The kagome 



antiferromagnet at low fields is described by the copla- 
nar triatic order parameter (H) with one of the spin triad 
vectors being antiparallel to the field direction, Fig. [l| 



H/Jn 




FIG. 1. Schematic phase diagram of the classical kagome 
antiferromagnet. The top inset shows the clapping mode for 
quasicollinear states. The lower inset shows a unit cell of the 
vo x v3 quasicollinear state: filled and empty circles denote 
the b- and the a-type spins, respectively. 

Generally, classical ground states at an arbitrary field 
< H < _ff sa t are all noncollinear. Collinear configura- 
tions appear among the ground state manifold only at 
special rational values of the applied field. The up-up- 
down (uud) structure has |Sa| = 1 and, hence, is stable 
only at H c = 2 J = ^H sat - The problem of calculating 
the total number of the uud states for the kagome lattice 
at H = H c can be mapped onto the problem of dimer 
coverings of a dual hexagonal lattice, which is solved ex- 
actly lijM . The number of uud states scales as V N , with 
V = 1.1137... A special role of such collinear spin config- 
urations in magnetic field for various frustrated models 
has been recently emphasized in M. Standard harmonic 
analysis for the kagome antiferromagnet at H = g-ff sa t 
yields three two-fold degenerate excitation branches for 
an arbitrary uud state: 



3 ± ^3(1 + 2 7k ) 



(5) 



with 7k = ^(cosfer + 2 cos \k x cos ^k y ). Therefore, the 
collinear states have _/V 4 = |A r soft quartic modes, twice 



more than any coplanar state. This fact has a simple 
geometrical origin. Local soft modes correspond to al- 
ternate tilting of spins around elementary hexagons. For 
collinear states such distortions have two polarizations in 
two directions perpendicular to the field, while coplanar 
states have soft modes only in the polarization transverse 
to the spin plane. 

Thermal fluctuations reduce the free energy of the uud 
states compared to other noncollinear classical ground 
states at H — ^H sat . Further lifting of degeneracy within 
the discrete subset of the uud states does not occur due 
to the same mechanism as for coplanar states at H = 
P,p|Jl4| . Every collinear state has special lines of weath- 
ervane defects, which contain alternating sequences of 



u-d-u-d spins. Flipping all spins along such a line by 
180° produces another collinear state at the free energy 
cost AF ~ T. A macroscopic temperature independent 
number of such defects is generated at low T's leaving 
all uud states accessible for the magnet and creating a 
collinear spin liquid. Enhanced number of soft modes is 
also reflected in the specific heat C: out of 2N degrees of 
freedom for N classical spins, every quadratic mode con- 
tributes \^b to the specific heat, while a quartic mode 
does only \k B §]. Hence, C = [l- N 4 /(4N)]k B per one 
spin and we expect C = Ifee for the collinear spin liquid 
compared to C = j^ks for coplanar states. Specific heat 
measurement in numerical simulations provides a useful 
way to distinguish different spin states. 

The free energy for the collinear spin liquid is lower by 
~ NT\n(J/T) compared to any noncollinear state. But 
a small deviation of field away from H = \H sat produces 
only a little change in the classical energy The collinear 
spin liquid is, therefore, stabilized by anharmonic inter- 
actions to a finite field range in a way similar to sta- 
bilization of the incompressible quantum fluid states in 
the Fractional Quantum Hall Effect. The collinear spin 
liquid phase yields a weak magnetization plateau at 1/3 
of the saturated magnetization: spins in the uud states 
cannot tilt towards the field and a weak variation of the 
magnetization for the collinear spin liquid is produced 
by thermal excitations alone. The above arguments also 
suggest the first-order transition from the collinear spin 
liquid state to the low- and the high-field phases at low 
temperatures. 

In the field range |-ff S at < H < if sat the ground state 
configurations are again all noncollinear. There appear, 
however, new quasicollinear states with two coinciding 
spins of the basis triad e^ = e c , which have zero chi- 
rality on every plaquette (Fig. |lj). These configurations 
belong to the subset of coplanar states. Hence, a qua- 
sicollinear state has all the soft excitations of coplanar 
states: ^ N quartic modes in the polarization transverse 
to the spin plane. In addition, they acquire extra in-plane 
soft modes. For a single plaquette, such a clapping-type 
mode corresponds to simultaneous tilting of the two par- 
allel spins in opposite directions inside the spin plane 
(Fig. nl). To preserve a weak anharmonic energy change 
the clapping deformation has to be extended on the lat- 
tice along a loop which contains only the 6-type spins. 
The shortest such loop is a perimeter of a hexagon, while 
number of appropriate hexagons is maximized if a quasi- 
collinear state is tiled in the \/3 x \/3 structure (Fig. nf). 
Thus, the \/3 x V3 quasicollinear state has the largest 
number of soft modes N4, = |iV at ^H sat < H < H sat . 

Nevertheless, the true breaking of the translational 
symmetry does not take place. The disordering mech- 
anism is again related to the presence of loops of alter- 
nating spin orientations or weathervane defects M . Loop 
flipping a-b-a-b — ► b-a-b-a at zero energy cost contributes 
~ T ln( J/T) to the thermodynamic potential because of 



a destruction of a few soft modes. Such losses are out- 
weighed at small concentrations of defects by the entropy 
gain. At finite temperatures there will be a macroscopic 
number of defects, which destroy the long-range v3 x v3 
translational symmetry and create instead a phase with 
a uniaxial triatic order parameter: 

S a ^ = rfnPrF - \ (n a 5( 3l + n% a + nH a p) , (6) 

where indices run over x,y and n = e a — (e a ■ z). 
Higher free energy costs make the number of defects in 
the uniaxial triatic state to be temperature dependent: 
Ndcf ~ TN with diverging magnetic correlation length 
£ ~ 1/T. Lacking a conventional magnetic order the 
high-field triatic phase asymptotically approaches the or- 
dered v3 x v3 quasicollinear structure. This behavior 
is reflected in the temperature dependence of the specific 
heat which should tend to C — |fes at low T. 




H/J 

FIG. 2. MC results for the classical kagome antiferromag- 
net. Top panel: susceptibility vs. magnetic field for several 
temperatures. The inset shows the magnetization curve at 
T = 0.01J. Bottom panel: specific heat vs. magnetic field 
for T = 0.002 J (full circles) and T = 0.01J (open squares). 
The inset shows the finite-size extrapolation for the squared 
triatic order parameter (g). 

At the saturation field H sa ,t = 6 J the ground state be- 
comes nondegenerate and corresponds to parallel align- 
ment of spins. The harmonic excitations in this collinear 
uuu phase are given by the same expressions ([3|) as for the 
uud states. There are N4 = |JV soft modes at this mag- 
netic field and the specific heat again reaches the value 
C — fkg. Soft modes also produce a universal magnetic 
behavior of the frustrated spin system at the saturation. 



For H > H sat all three branches (g) acquire an additional 
field dependent shift Su> — (H — H sat ). Taking into ac- 
count only the most singular contribution of the first flat 
branch, we derive that at H > -ff sa t the classical partition 
function scales as 



Z(H, T) = T l ' A f{u) , u = (H- H sat )/VT 



(7) 



The susceptibility x = dM/dH determined from M) is 
X = [/"(«)//(«)] - [f'{u)/f{u)f. Atu = 0(H = ff sat ) 
the susceptibility becomes temperature independent and 
all the curves x(-ff, T) cross at one point, see Fig. ^. 

So far we considered only the short wave-length degrees 
of freedom and completely disregarded long wave-length 
excitations, which determine the nature of phase tran- 
sitions in two dimensions. The non-Abelian topological 
defects discussed previously for the triatic order parame- 
ter H are suppressed by the field. The two triatic phases 
break only U(l) symmetry for rotations about the field 
direction. Berezinskii-Kosterlitz-Thouless type transi- 
tions separate them from the high temperature paramag- 
netic phase, see Fig. [|. In contrast, the first-order tran- 
sition to the collinear spin-liquid state survives and cor- 
responds to a gas-liquid type transition on the boundary 
with the paramagnetic phase. The collinear uuu phase 
also has an enhanced entropy contribution and is stabi- 
lized below Hsat- The corresponding line of first-order 
transitions ends at a critical point as sketched in Fig. [I]. 

We have also performed Monte Carlo (MC) simula- 
tions for the model (pj). The standard Metropolis algo- 
rithm was used with up to 10 MC steps per every point. 
Simulations were done for periodic 3L 2 -site lattices with 
L = 6, 12, 18; the data presented in Fig. are for a 
972-spin cluster. The magnetization curve at T = is 
a straight line with a slope x — dM/dH = \- At low 
temperatures x{H) develops a dip near H = \H sat corre- 
sponding to a weak plateau at i-M sa t. Two peaks, which 
surround the dip, indicate first-order transitions to the 
plateau phase. The peaks become rounded and, then, 
completely disappear at higher temperatures. We esti- 
mated T* ~ 0.025J as the highest temperature, where 
the collinear spin-liquid still exists (Fig. |l|). Above this 
temperature the dip corresponds to a smooth crossover. 
The field dependence of the specific heat also follows the 
above predictions. At T = 0.002 J, the specific heat starts 
near y^fes, corresponding to the coplanar triatic state, 
at low fields. As the field approaches \H sa , t , C{H) goes 
down to §/cb, corresponding to the collinear spin liq- 
uid, and, then, recovers back to §£;_b, corresponding to 
the uniaxial triatic state. The specific heat drops again 
near the saturation with C|# sat = |fcs, while the min- 
imum value is reached at somewhat lower fields. Such 
a behavior signals stabilization of the uuu phase beyond 
its classical boundary: bare negative modes below -ff sa t 
are renormalized into soft modes, which further decrease 
C(H). The specific heat also exhibits a peak at H ss 4.8 J 



(T = 0.002J) and H « 3J (T = 0.01J) with a signifi- 
cant finite-size dependence. We attribute this peak to a 
phase transition from the collinear uuu state into a non- 
collinear uniaxial triatic state (g) . To check this we have 
calculated field dependence of the squared triatic order 
parameter (0) at T = 0.002 J for three cluster sizes and 
extrapolated it to the thermodynamic limit. Results are 
shown in the inset in Fig. 0. The uniaxial triatic phase 
disappears at about the same field as the position of the 
peak in the specific heat. The present data do not resolve 
clearly the order of the phase transition, which must be 
determined in a detailed numerical study. 

We have shown that the geometrical approach based 
on soft mode counting is a powerful tool for investiga- 
tion of the field behavior of a classical kagome lattice 
antifcrromagnet. Our preliminary analysis shows that 
a Heisenberg antifcrromagnet on the related garnet lat- 
tice has the same type of the phase diagram, Fig. [y. In 
particular, the field range of an asymptotically ordered 
quasicollinear state coincides roughly with a region for a 
field-induced ordering observed in GdaGasO^ ||. 
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